On the gravitational and electrodynamical stability of massive nuclear density cores 
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We present a unified treatment of nuclear density cores recovering the classic results for neutral 
atoms with heavy nuclei having a mass number A ~ 10 2 -10 6 and extrapolating these results to 
massive nuclear density cores with A « (mpi anc k/m n ) 3 ~ 10 57 . The treatment consists of solving the 
relativistic Thomas-Fermi equation describing a system of N n neutrons, N p protons and N e electrons 
in beta decay equilibrium. The N p protons are distributed at a constant density within a spherical 

/ \ 1/2 3 

core of radius R c . A new island of stability is found for A > A R = 0.039 ( J ( m "° T1 ° ck ) 



The 



Coulomb repulsion, screened by relativistic electrons, is balanced by the gravitational self-interaction 
of the core. In analogy to heavy nuclei they present, near their surface, an overcritical electric field. 
The relation between A and N p is generalized to an arbitrary value of the mass number, and the 
phenomenological relations for A < 1.5 • 10 2 are obtained as a limiting case. 
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The study of neutral atoms with heavy nuclei of mass 
number A ~ 10 2 — 10 6 is a classic problem of theoretical 
physics 0, Special attention has been given to the 
study of their stability against Coulomb repulsion [3J] and 
to the existence of electric fields larger than the critical 
value 

E c = m 2 e c 3 /eh, (1) 

near their surfaces We have extrapolated these 

results by numerical integration to the case of mas- 
sive nuclear density cores: an overall neutral system 
of neutrons, protons and electrons, in beta decay equi- 
librium, at nuclear density and having mass numbers 
A ~ (mpi anc k/r7i n ) 3 where m n is the neutron mass and 

mpia n ck-(VG) 1/2 1- 

In this letter we present a unified treatment of heavy 
nuclei and of massive nuclear density cores made pos- 
sible by an explicit solitonic solution of the relativistic 
Thomas- Fermi equation which governs these phenomena. 
We confirm the existence of overcritical electric fields 
near the surface of massive nuclear density cores. The 
maximum value is given by E max — 0.95y/am 2 c 3 / eh 
where as usual a — e 2 /(he) and is the pion mass. In 
contrast with the case of heavy nuclei we prove the exis- 
tence of a new island of stability for mass numbers obey- 
ing the condition A > A R = 0.039 (^f) V2 ( m ^ a ; ck ) 3 - 
In this case the equilibrium against Coulomb repulsion 
is not due to the surface tension of strong interactions. 
It originates from the over-whelming effects of the self- 
gravitational interaction of the massive dense core. Fi- 
nally we obtain a generalized relation between the mass 
number A and atomic number N p which encompasses 
previous phenomenological expressions. 

It is well known 0] that the stability of nuclei in neutral 
atoms is guaranteed by the surface tension £ s w 17.5 • 



A 2 / 3 MeV, created by the nuclear forces balancing the 
Coulomb repulsion £ cm w (3/5)e 2 N 2 /R c , where R c — 
TqA 1 / 3 is the nuclear radius and r 
0.8b(h/m v c). If one assumes 



1.2 • 10~ 13 cm 



v 2 ' 



(2) 



and a constant proton density n p = N p hq/A 
0.25(m 7r c/?i) 3 , where no = A/(4ttR 3 /3) is the ordinary 
nuclear density, £ s > £ em gives the stability condition 
A < A* — 10~ 1 m„/(m 7r a) [§]. If one assumes (see (lf3 |) 
a more accurate phenomenological expression relating N p 
and A 
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A 200^/3 



(3) 



the stability condition is somewhat cumbersome and cor- 
responds to A < A* w 1.5 • 10 2 . 

A novel situation occurs when super-heavy nuclei (A > 
A ~ 10 4 ) are examined @, ■ The distribution of elec- 
trons penetrates inside the nucleus: a much smaller effec- 
tive net charge of the nucleus occurs due to the screening 
of relativistic electrons @, [13, 11 1 ■ A treatment based on 
the relativistic Thomas-Fermi model has been developed 
in order to describe the penetration of the electrons and 
their effective screening of the positive nuclear charge. 
In particular, by assuming Eq. ((2|), in a series of classic 
papers Greiner et al. [H, 0, 01 and Popov et al. [3, H, 0] 
were able to solve the non-linear Thomas-Fermi equa- 
tion. It was demonstrated in that the effective posi- 
tive nuclear charge is confined to a small layer of thickness 
~ h/y/am^c. Correspondingly electric fields of strength 
much larger than the critical value given by Eq. ([T]) for 
vacuum polarization at the surface of the core are cre- 
ated. Under these conditions, however, the creation of 
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electron-positron pairs due to the vacuum polarization 
process does not occur because of the Pauli blocking by 
the degenerate electrons [12j. Although the screening of 
the positive charge of the core increases with increasing 
N p , the super-heavy nuclei nevertheless remain unstable 
against the Coulomb repulsion Q. 

Here we generalize the classic work of Greiner and 
Popov to the case A w {mp\ anc ^/m n ) 3 ~ 10 57 . We 
have also relaxed the condition expressed in Eqs. (J2J> they 
adopted, by explicitly computing the beta decay equilib- 
rium between neutrons, protons and electrons. A super- 
critical field still exists in a shell of thickness ~ Tij ' ^/am^c 
at the core surface, and a charged lepton-baryonic core 
is surrounded by an oppositely charged leptonic com- 
ponent. Such massive nuclear density cores are glob- 
ally neutral. We show that they are stable against the 
Coulomb repulsion of the proton component due to the 
stabilizing effects of the gravitational self-interactions. 

The analytic solution representing a core of degenerate 
neutrons, protons and electrons is obtained by assuming 
its radius to satisfy 



Rr 
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(4) 



where A is a parameter such that A « 1 (A < 1) corre- 
sponds to nuclear (supranuclear) densities when applied 
to ordinary nuclei. The overall Coulomb potential satis- 
fies the Poisson equation 



V 2 V(r) = -4ttc [n p (r) - n e (r)] , 



(5) 



with the boundary conditions V^(oo) = (due to the 
global charge neutrality of the system) and finiteness of 
V(0). The density n e (r) of the electrons of mass m e and 
charge — e is determined by the Fermi energy condition 
on their Fermi momentum Pf 

E F e = [(Pfc) 2 + my] 1 ' 2 - m e c 2 - eV{r) = , (6) 



which leads to 
(Pf) 3 



n e (r) 



1 



37T 2 h 3 3tt 2 ^c 3 



, [e 2 V 2 (r) + 2m e c 2 eV(r)] V2 . 



(7) 



By introducing the dimcnsionless quantities x — 
r /[h/m-n-c], x c = R c /[h/m„c] and x/ r = eV(r)/ch, the 
relativistic Thomas-Fermi equation takes the form 



1 d 2 X (x) 
3x dx 2 



X 2 (x) 



3/2 



which in turn is related to the proton and electron densi- 
ties by Eqs. ([5]), Q and ((SJ). These equations have been 
integrated numerically ||. 

In the ultrarelativistic limit, the relativistic Thomas- 
Fermi equation admits an analytic solution. Introducing 
the new function (j> defined by 



(8) 



where x(0) = 0, x(°°) — 0- The neutron density n n (r) is 
determined by the Fermi energy condition on their Fermi 
momentum Pi imposed by beta decay equilibrium 

El = [{Plc) 2 + m 2 n cV 2 - mn c 2 

= [(Pic) 2 +m 2 c 4 } 1 / 2 -m p c 2 +eV(r), (9) 
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and the new variables x = (12/n) \/aA 1 x, £ — x—x c 
where x c = (12/ir) 1 ^ ^fa~A~ 1 x c , then Eq. © becomes 



d 2 m 



(10) 



where <j)(£) = </>(£ + x c ). The boundary conditions on 
<j> are: (/)(£) — > 1 as £ — > — x c -C (at the massive nu- 
clear density core center) and </>(£) — > as £ — > oo. The 
function 4> and its first derivative 4>' must be continuous 
at the surface £ = of the massive nuclear density core. 
Eq. (TIT)]) admits an exact solution 



(11) 



f l-3[l + 2- 1 / 2 sinh(a-V30] t < 0, 

= { J2_ " . „ 

(£ + &)' 



£>0, 



where the integration constants a and b have the values 
a = arcsinh(ll % /2) = 3.439, b = (4/3)V2 = 1.886. We 
can next evaluate the Coulomb potential energy function 



'Qtt\ 1/3 1 
eV(0 = [—) £"^0(0, 



and by differentiation, the electric field 
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(13) 



Details are given in Figs. [I] and [2] 

Next we can estimate two crucial quantities: the 
Coulomb potential at the center of the configuration and 
the electric field at the surface of the core 



0.95y/a- 
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(14) 



E c . (15) 



Remarkably these two quantities are functions only of the 
pion mass m^, the density parameter A and of course the 
fine constant structure a. Their formulas apply over the 
entire range from superheavy nuclei with N p ~ 10 all 
the way up to massive cores with N p rs (mpi anc k/?Tin) 3 . 

Using the solution (fTTj) , we have obtained a new gen- 
eralized relation between A and N p for any value of A 
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S, Mass number A 



FIG. 1: The electron Coulomb potential energy — eV, in units 
of pion mass m-n is plotted as a function of the radial coordi- 
nate £ = x — x c , for selected values of the density parameter 
A. 
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FIG. 2: The electric field is plotted in units of the critical field 
E c as a function of the radial coordinate £ for A=2, showing 
a sharp peak at the core radius. 



which agrees remarkably well with the phenomenological 
relations given by Eqs. ^ and © in the limit A < A*, 
as clearly shown in Fig. [3] 

Having established the validity of the above equations 
both for superheavy nuclei and for massive nuclear den- 
sity cores, we now outline some fundamental differences 
between these two systems. The charge-to-mass ratio of 
the effective charge Q at the core surface to the core mass 
M is given by 



Q 



GM 



Gm n A 



WlPlanck 
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Np 
A 



(16) 



FIG. 3: Our A-N p relation at nuclear density (solid line) 
obtained from first principles compared with the phenomeno- 
logical expressions given by Eqs. ([2]) (dashed line) and (|3]) 
(dotted line). The asymptotic value, for A — > (mpi anc k/m„) 3 , 



is N p « A/ 38 (see J3) 



mass ratio for the nucleus is 

Q 1 TOpianck 



GM 20 



10 



18 



(17) 



There is obviously no role in the stabilization of these 
nuclei by the gravitational interactions. 

Instead for massive nuclear density cores where N p « 
(wpianck/'Tin) 3 , the ratio Q/VGM given by Eq. (|T6]) is 
simply 

7m~i^ (18) 

which is approximatively 1/38 (see Fig. [3]). It is well- 
known that the condition that the charge-to-mass-ratio 
(|18p be smaller than 1 is a necessary one for the equi- 
librium of self-gravitating mass-charge system both in 
Newtonian and general relativity (see, e.g., [l]| ). Thus 
massive nuclear density cores are gravitationally stable, 
globally neutral and bound. It is therefore possible to 
formulate for them a consistent stable model in terms of 
gravitational, strong, electro mag netic and weak interac- 
tions and quantum statistics [lq . 

We can see the gravitational stability of massive nu- 
clear density cores from a different point of view. The 
maximum Coulomb energy per proton is given by Eq. I|14p 
where the potential is evaluated at the center of the core. 
The gravitational potential energy per proton (of mass 
m p ) in the field of a massive nuclear density core with 
A w (m P i anck /m„) 3 , is given by 



£ g = 



c Mm p 

Rr. 



1 ^Planck rn-^c 2 



N, 
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For superheavy nuclei with N p w 10 3 , the charge-to- 



1/3 



(19) 
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Since A/N p ~ 38 (see Fig. [3]), independently of A value, 
the gravitational energy is larger in magnitude than and 
opposite in sign to the Coulomb potential energy per pro- 
ton of Eq. (TT4"]) so the system is gravitationally stable. 

There is yet a third more accurate derivation of the 
gravitational stability based on the analytic solution of 
the Thomas- Fermi equation Eq. (flQ| . The Coulomb en- 
ergy £ em is @ 



E 2 
8tt 
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R 2 c (eV(0)) 3 
{Zna) 1 / 2 
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(20) 



which is mainly distributed within a thin shell of 
width 5R C ~ hA/(^/am w c) and proton number 5N P = 
npAnR^SRc at the surface. To ensure the stability of the 
system, the attractive gravitational energy of the thin 
proton shell 



_ G Mm p 6N p 



Rc 
t GA*/ 3 
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(21) 



has to be larger than the repulsive Coulomb energy (|2"0]l. 
For small A the gravitational energy is always negligible. 
However, since the gravitational energy increases propor- 
tionally to A 4 / 3 and the Coulomb energy proportionally 
to A 2 / 3 , such a crossing necessarily exists. We obtain the 
crossing at 

1/2 
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This establishes a lower limit for the mass number Ar 
necessary for the existence of an island of stability for 
massive nuclear density cores. 

Thus the arguments often quoted, concerning limits on 
the electric fields of an astrophysical system based on a 
free test particle approximation given by equations like 
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(23) 
(24) 



appear to be inapplicable for A ~ {mpi anc k/m n ) 3 , when 
the collective effects of the quantum statistics are present 
and properly taken into account through the relativistic 
Thomas-Fermi model. Eqs. (|23|) and (f24|) have to be 
replaced by Eqs. JTSJ) and (fl"8|). 
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Having established the role of gravity in stabilizing the 
Coulomb interaction of the massive nuclear density core, 
we outline the importance of the strong interactions in 
determining its surface. We find for the neutron pressure 
at the surface: 
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(27) 



and for the surface tension, as extrapolated from nuclear 
scattering experiments, 



P, 



O.m ( A\ 2 ' 3 1 

4^ J {h/m^cf \N P ) A 2 ' { ' 



We then obtain 



J§i = 0.39 • A 3 ( ^ ) = 0.24 



Pnucl 



(29) 



dust 



where p nuc i = m n nQ. The relative importance of the 
nuclear pressure and nuclear tension is a very sensitive 
function of the density at the surface p su rf i which in the 
present case is also constant throughout the configura- 
tion. This ratio will in turn determine the extent of a 
"crust" surrounding the core, with leptonic pressure and 
with density due to nuclei. 

In conclusion, we have obtained the following results. 

a) We have found a new island of stability in addi- 
tion to the known one when the Coulomb repulsion of 
heavy nuclei is balanced by the surface tension of strong 
interactions Q . The new island occurs for much larger 
mass numbers: for A > Ar as given by Eq. (|22p . The 
Coulomb repulsion is now balanced by the gravitational 
forces and the role of strong interactions is only relevant 
in balancing the neutron pressure and determining the 
size of the core. 

b) The systematic use of the relativistic Thomas- Fermi 
equation and the enforcement of beta equilibrium has al- 
lowed to obtain, from first principle, a relation between 
the mass number A and the proton number N p both in 
the mass range of stable nuclei and in the massive nu- 
clear density cores for A w (mpi anc k/m n ) 3 ~ 10 57 . For 
A < A* this new relation approaches asymptotically the 
phenomenological relations expressed by Eqs. @ and 
©. For A > A R it leads to N p /A - 1/38. 

c) Although the configurations are globally neutral, 
electric fields, whose order of magnitude is larger than 
the critical field given by Eq. ([1]), are confirmed to exist 
close to the surface of the core. The charge to mass ratio 
at the surface, as well as the values of the electric field 
are of the order of 10 14 larger than the ones expected 
from simpler dust approximations. 

We are currently relaxing the constant proton density 
condition adopted in this article taking also into account 
the general relativistic effects by a generalized Tolmann- 
Oppenheimer-Volkoff equation 15j in order to describe 
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the gravitational interactions of such cores . The ap- 
plications of these results to the physics of neutron stars 
appear to be particularly promising H| • We are also 
exploring the effects of the electromagnetic structure of 
these massive nuclear density cores during the process of 
gravitational collapse to a Black Hole. 
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